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Abstract 

We investigate the behavior of the zero-temperature quantum non-linear sigma 
model in d dimensions in the presence of a damping term of the form f(uS) ~ M", 
with 1 < a < 2. We find two fixed points: a spin-wave fixed point FP1 showing a 
dynamic scaling exponent z = 1 and a dissipative fixed point FP2 with z > 1. In the 
framework of the e-expansion it is seen that there is a range of values ct*(d) < a < 2 
where the point FP1 is stable with respect to FP2, so that the system realizes a 
2 = 1 quantum critical behavior even in the presence of a dissipative term. However, 
reasonable arguments suggest that in d = 2 this range is very narrow. In the broken 
symmetry phase we discuss a phenomenological scaling approach, treating damping 
as a perturbation of the ordered ground state. The relation of these results with the 
pseudogap effect observed in underdoped layered cuprates is discussed. 

Key words: renormalization group, non linear sigma model, quantum critical 
point, cuprates. 
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1 Introduction 



The presence of an antiferromagnetic (AF) phase rapidly replaced by an 
anomalous metallic phase upon doping is one of the prominent features of 
the phase diagram of high temperature superconducting cuprates [1,2]. In- 
deed it was early suggested [3] that the proximity to an insulating magnetic 
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phase together with the nearly twodimensional structure of these materials 
could be responsible for the superconducting and anomalous normal-state (i.e. 
non-Fermi liquid) properties of these systems via the creation of a quantum- 
disordered spin liquid (resonating- valence-bond state). Since then a great deal 
of attention has been devoted to the fascinating interplay between magnetism 
and charge degrees of freedom, leading to various theoretical proposals and to 
different scenarios for the metal-insulator transition at low doping. 

One scenario is based on the tendency of a magnetically ordered phase to seg- 
regate the additional holes due to doping [4] . The tendency to phase separation 
is then frustrated by the Coulombic repusion between the segregated carriers 
[5], thus leading to formation of hole-rich domain walls separating AF domains 
in the form of stripe textures. Within this scenario, pseudogaps naturally arise 
in the underdoped phase of the materials as a consequence of stripe fluctua- 
tions and local pair formation [6]. Obviously this will substantially affect the 
transition to the AF insulating phase at low doping [7,8]. 

The transition between a magnetically ordered phase and a metallic state can 
also be strongly affected by extrinsic ingredients like the disorder induced by 
dopant ions (Sr or Zn in the case of La2_ x _ y Sr x Zn y Cu0 4 ). This may lead to the 
formation of local random magnetic moments giving rise to spin-glass ordering 
between the AF insulating and the paramagnetic metallic (PM) phases [9- 
12,1]. 

The stripe and the spin-glass scenarios clearly illustrate the complicated and 
rich nature that the AF-PM transition may take. However, other proposals are 
more directly related to magnetism and to the contiguity between the AF and 
the metallic phase. These proposals start from the seminal work of Ref. [13] 
showing that the materials with half-filled Cu0 2 planes are suitably described 
in terms of a twodimensional quantum Heisenberg AF model with very small 
interplanar coupling. Following work [14-18] suggested that, although the 2D 
model displays long-range order at T = 0, upon doping the added charges 
enhance the quantum spin fluctuations thus driving the system into a disor- 
dered state even at zero temperature (quantum disordered phase) [14]. The 
occurrence of an AF quantum critical point (QCP) provides a natural frame- 
work to intepret the scaling behavior at low temperatures and low frequen- 
cies of the of the g-integrated magnetic susceptibility x{°°) = I ^qxiq,^)- 
The susceptibility x{ u ) can be experimentally determined via magnetic res- 
onance measurements of the nuclear copper spin-echo decay rate, T 2 q and 
spin-lattice relaxation rate, Tf 1 [19,2]. It is seen that in a range of tempera- 
tures T* < T < T cr , depending on doping, T 1 T/T 2 g = const [20]. Assuming the 
scaling relations [14], x(Q,u) ~ t 2 ~ v x(\q ~ Qaf\^u/T) and f ~ T" 1 / 2 , with 
z a dynamic critical exponent, r] an anomalous dimension, Qaf the antiferro- 
magnetic ordering wave-vector, £ the correlation length, x a scaling function, 
one gets T{I '/T 2G ~ T 1 ^ 1 / 2 [17,18]. The constancy of this ratio is thus inter- 
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preted as the signature of a critical behavior corresponding to a z — 1 value 
of the dynamic index [17,18,21-23]. The strong critical spin fluctuations oc- 
curring at the AF-QCP have also been claimed to provide a possible source 
of pairing and normal-state anomalies. 

Despite the apparent simplicity of a scenario involving a direct transition 
between an insulating AF phase and a PM phase, such a transition is far from 
being trivial and is still an open problem. The main difficulty is in the lack of 
a microscopic model being able to smoothly interpolate between an insulating 
AF ordered phase with gapped charge excitation and a PM phase. On the 
one hand, the insulating AF phase is suitably represented by the Heisenberg 
model and, in particular, the long- wavelength and low frequency behavior of 
spin fluctuations in the Cu0 2 planes can be well described by the quantum 
non-linear sigma model [13], which at zero temperature and in two spatial 
dimensions is characterized by the action 

S = 1 J d 2 k (k 2 + ^ ■ <t>_ k ^ (1) 

Here 4> is a three-component vector field subject to the real-space constraint 
<j) 2 X T = 1 which describes the local staggered magnetization, and g is a coupling 
constant. 

It was first argued in [14] that this description could be carried over to the case 
of the lightly doped compounds. In fact, one can think that in this case the 
4> field continues giving an effective description of the electron spins localized 
on the Cu sites of the lattice, while the presence of itinerating holes provides 
a finite renormalization of the coupling constant. The theory described by the 
action (1), besides having the correct symmetries, realizes a zero-temperature 
transition for a given critical value g c of the coupling (see [13,24] for the 
quantum version and [25-27] for the classical theory) thus providing a good 
model for the physical picture of the AF-QCP [15]. From this point of view 
the z — 1 value of the dynamic scaling exponent is a simple consequence of 
the Lorentz invariance of (1). 

However, it is clear that the non-linear a model does not properly describe 
the low-energy charge excitations, which are characteristic of a metallic state 
possibly providing a damping mechanism for spin- waves (see, e.g., case B of 
Ref. [28]). Therefore the model in Eq. 1 provides a good description of a 
transition between an insulating (charge-gapped) AF (AFI) and an insulating 
paramagnet (PI)., but it appears to be inappropriate when spin- waves can 
decay in low-energy particle-hole excitations. 

On the other hand, attempts have also been made to provide a microscopic 
basis for the transition to a PM phase. In this regard two distint cases A and 
B depending on the shape of the Fermi surface in the metal were investigated 
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in Ref. [28] within a spin-density-wave approach. In case A the AF wavevector 
Qaf = (*7r/a, 7r/a) (a is the lattice spacing) was not able to connect points of 
the Fermi surface, so that momentum and energy conservation did not allow 
the decay of spin-waves in particle-hole pairs. In this case the transition was 
naturally characterized by a dynamic critical index z — 1. In case B, instead, 
there were ("hot") points on the Fermi surface connected by Qaf so that 
low-energy particle-hole pairs could be created by the decay of low-energy 
spin-waves. However, in both cases, A and B, the phase with AF long-range 
order is metallic and therefore the model does not suitably describe a transition 
between an AFI and a PM The case A above shares some common features 
with a mixed spin-fermion model without a three-body (i.e. Kondo-like) direct 
coupling between the local moments and the fermion spins [29]. This model 
has been shown [30] to display a T = transition between a metallic (Fermi 
liquid) phase with spin commensurate long-range order and a metallic (also 
Fermi liquid) quantum disordered phase. Owing to the absence of the three 
body coupling responsible for the direct spin-wave decay into particle-hole 
pairs, also the QCP of this model is characterized by z — 1. 

To the best of our knowledge no simple microscopic modelization is presently 
available to describe a direct AFI-PM transition. Of course the possibilities 
remain open that other extrinsic mechanisms (stripes, disorder, ...) provide 
intermediate steps between the AFI and the PM phases. Alternatively, one 
may also envisage that an AFI to a PI phase transition occurs first at lower 
doping, followed at slightly larger doping by a PI-PM transition between two 
spin-disordered phases. 

Together with the basic difficulty outlined above, the additional question arises 
concerning the observed z = 1 critical behavior in the underdoped phase of 
the cuprates and the way this can be reconciled with the metallic character of 
this phase. In particular it is believed that the Fermi surface of the cuprates 
in the metallic phase contains "hot" points connected by the AF wavevector. 
Therefore it should be possible for the spin-waves of arbitrarily low energy 
to decay into particle-hole pairs and get overdamped. This would lead to a 
relaxational behavior contrasting with the z = 1 (i.e. propagating) critical 
behavior of the spin-waves. Moreover, by accepting a z = 2 critical behavior 
in the metallic phase, it remains to clarify how this may turn into the natural 
propagating behavior, which characterizes the spin excitations in the ordered 
side of the QCP. This latter difficulty was also encountered and stressed by 
the authors of Ref. [28] for their case B: within their mean-field description 
a z = 2 behavior was obtained both in the disordered and the ordered AF 
(metallic) phases. This suggested that the simple formulation B of the mean- 



2 However, owing to the specific form of the spin-wave damping for a Fermi surface 
of type A, the non-linear a model would suitably describe the system also in this 
case despite the metallic character of both the ordered and disordered phases. 
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field SDW model was only applicable on the disordered side of (and not too 
close to) the AF-QCP. 

In the present paper we will not attempt to find a microscopic foundation to 
a direct AFI-PM transition, but we will instead start from a semiphenomeno- 
logical point of view by introducing a field-theoretical generalization of the 
non-linear a model action (1) to investigate a) the condition to be fullfilled by 
a simple (fc-independent) damping term in order to maintain a z — 1 critical 
behavior; b) how in the z > 1 case an additional crossover energy must behave 
in the ordered phase to provide a scale separating damped spin excitations at 
high energy and propagating spin excitations at low energy. 

The resulting theory will depend on the precise form of the effective damping 
term. We first discuss the simpler case of the metallic phase on the g > g c 
side of the QCP. Within a microscopic model of spins coupled to free itinerant 
holes, integrating out the holes degrees of freedom in the random phase ap- 
proximation gives a damping term of the form f(u) = j\uj\ (see e.g. [31,28]). 
However, the choice of the adequate low frequency form of f(u>) should also 
take in consideration the fact that in underdoped cuprates a significant loss 
of low-frequency spectral weight is observed at low temperatures in the distri- 
bution of quasi-particles, a phenomenon often referred to as "pseudogap" (see 
e.g. [32]). The origin of this effect is a debated issue, and has been alterna- 
tively interpreted: as the result of spin-fermion scattering [33], as due to the 
intrinsecally non- Fermi liquid nature of the system (Luttinger liquid) [34], as 
a signal of the formation of preformed pairs for T > T c [32,35], as due to stripe 
charge fluctuations [6]. In all cases, a precise computation of these effects im- 
plies solving a difficult self-consistency problem going well beyond the RPA 
approach. We shall try here a rough, phenomenological approach, assuming 
that the damping coefficient 7 can acquire at low frequencies an cu-dependence 
~ \ou\ x , thus giving rise to an effective damping term of the form f(u>) ~ \&\ a , 
with a = l + x>lQ. 

Once chosen a damping term f(w) we perform a renormalization group (RG) 
analysis of the resulting theory in the proximity of the QCP using the mo- 
mentum-shell method and e-expansion. We observe the existence of two fixed 
points, which we shall denote FP1 and FP2, showing respectively a z = 1, 
undamped spin- wave behavior and a z > 1, dissipative behavior. The point 
FP1 corresponds to the zero-temperature critical point of a pure quantum 
non-linear sigma model without any damping term, but it becomes unstable 
in a wide region of the parameters with respect to the dissipative fixed point 
FP2. 

There exists a range of values a* < a < 2 for which one obtains a z = 1 

3 The a = 1 case corresponds to the case considered in Ref. [28] within a large-N 
framework. 
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dynamic scaling exponent in the proximity of the stable quantum critical 
point, even in the presence of a damping term. In other words, there are "soft" 
damping terms which do not destroy the z = 1, spin- wave behavior. Here a* is 
the dimension dependent lower bound for the exponents in the soft damping 
term, which we determine to be a* = 2 — 77(0;*) = 2 — 771, where 77(a) is the 
anomalous field dimension and 771 is its value at the FP1. It is worth noting that 
this result b resemblance to what is found in classical models with long- 
range interactions [36,37]. In particular the scaling exponent of the damping 
stems from the non-analytic form of f(ui), which is not reconstructed by the 
RG corrections (for a discussion of this point in classical models see Ref. [38]). 
Therefore 7 does not acquire singular corrections but for the anomalous field 
dimensions. 

Within the framework of an expansion in e = d + z^°> — 2, where d is the 
spatial dimension of the system and z<® is the bare value of the dynamic 
scaling exponent in the vicinity of the stable fixed point, one would get 
a* = 2 — (ae 2 )/(l + a) = 2 — ei, where the 1 or 2 subscripts in e refer to 
FP1 or FP2 respectively. However, the extrapolation of results valid for e ~ 
to the physical region d = 2 should be done with great care. In particular, 
using the numerical results for the critical exponents of the 0(3) model in 2+1 
dimensions, it appears that a* is very near to two and the range of "soft" a val- 
ues should be a very narrow one. This means that within the here-considered 
model the appearance of a z = 1 critical behavior in the cuprates cannot be 
accounted for in the effective framework without assuming a value a ~ 2, i.e., 
an almost linear decay of the quasi-particle spectral weight for low frequencies. 
This spectral weight suppression is a rather severe condition to be obtained 
from the magnetic scattering itself, as also suggested by a direct perturbative 
calculation within a spin-fermion model [39] , where density of states and ver- 
tex corrections to the fermion bubble produced only a minor change of the 
power-law dependence of the damping term. It is worth noting that a substan- 
tial spectral weight reduction is however consistent with Anderson's idea that 
a smooth connection between a charge-gapped insulator and a Fermi-liquid 
metal is hardly conceivable. 

As far as the ordered phase is concerned, (point (b) above) it is necessary 
to take into account the presence of a non-zero staggered magnetization Nq, 
which is accompanied by the presence of two quasi-particle bands separated 
by an energy A. This suggests [28] to introduce a damping term having two 
distinct asymptotics: f(u) ~ \u\ a for \u\ ^> A and f(u) ~ ui b for |a>| <C A, 
with l<a<2<6<30 The difficult point here is that the scale A 



4 In principle the damping could also depend on q = Q AF +fc in the ordered (as well 
as in the disordered) phase, but in a first approximation we ignore this dependence 
mainly for simplicity reasons, but also because in the insulating phase, where a 
gap for the charge excitations providing spin- wave decay is present, it would be 
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separating damped and undamped spin excitations vanishes at the transition. 
This prevents using a conventional RG approach, and we will resort to simple 
scaling arguments and phenomenological assumptions. Our analysis does not 
attain the same degree of reliability of a traditional Landau-Ginzburg- Wilson 
theory and has a heuristic character. We shall assume that A has a power-law 
behavior, to be determined by matching the hydrodynamic behaviors on the 
two sides of the quantum transition. For a < a*, i.e. for z > 1, we find that 
(i) A should close at the "microscopic" scale as £ _z for £ — > oo, where £ is a 
correlation length, and that (ii) the spin- wave velocity c sw vanishes as to 
be consistent with the presence of damping. 

The paper is organized as follows: in Sec. 2 we derive the RG equations for 
the zero-temperature quantum non-linear sigma model with a damping term 
of generic form, following the scheme exposed in [13,40]. This part is rather 
technical and it can be skipped in a first reading without loosing the general 
meaning of the paper. 

In Sec. 3 we perform the fixed point analysis of the RG equations for a damp- 
ing term f{uj) = 7|c^| a . Here we show the existence of the two fixed points 
mentioned above, study their respective regions of stability and compute the 
relevant critical exponents, with particular regard to the value of the dynamic 
scaling exponent z. 

In Sec. 4 we discuss damping in the ordered phase, introducing a phenomeno- 
logical scale A vanishing at the quantum critical point with a critical exponent 
that we determine by matching the hydrodynamic behaviors on the two sides 
of the transition. We also compute the critical behavior of the spin-wave ve- 
locity c sw in the presence of damping. 

In Sec. 5 we present some conclusive remarks and discuss the connection be- 
tween our model and other physically interesting systems, like macroscopic 
quantum tunnelling and spin chains with long range interaction. 



2 Renormalization group equations 

Let us consider the Euclidean action of the quantum non-linear sigma model 
[13,24-26] at zero temperature in d spatial dimensions, with a damping term 
f{uj) of generic form and in the presence of a constant magnetic field h: 



S=- fdK + ^ 

2gJ \ c 2 




K 



(2) 



unnatural to expect a relevant q-dependence. 
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Here </> is an iV-component field satisfying the constraint <f> = 1 in real space, 
k = (k,uj), dn = (27r)~^ +1 ^ d d kduj. The magnetic field is oriented along the 
iV-th direction. The momentum and frequency integrations have been rescaled 
so that the adimensional spin velocity c has bare value equal to one and the 
momentum integration is cut off at A = 1. 

With our rescaling, in (2) the dimensionless coupling constant g is related 
to the bare (dimensional) spin- wave velocity Cq, the physical cut-off A and 
the bare spin-stiffness constant p by the relation g = hc A d ~ 1 /pQ [13]. In 
the large S (here S is the spin) limit p an d c are related to the coupling 
constant of the Heisenberg model J and the lattice spacing a by p = JS 2 a 2 ~ d 
and c = 2\fdJSa/h [24] with A ~ a .his also adimensional. 

Action (2) is formally equivalent to the Hamiltonian of an anisotropic corre- 
sponding statistical mechanical model in d + 1 dimensions. The critical point 
g — of the classical model is here reinterpreted as a zero-temperature quan- 
tum critical point, separating the ordered and the disordered ground state of 
the system, which are realized, respectively, for g < g* and g > g*. In order to 
study the long-distance, low-frequency behavior of the 4> field in the proximity 
of the quantum critical point we shall use the renormalization group and the 
e-expansion, following Refs. [13,40]. 

In the broken symmetry phase we let 7r = ((f) 1 , . . . , Ar ~ 1 ), a = <p N . Keeping 
only terms up to 0(g 2 ), action (2) takes the form 



+ 2# / d(4)/t i~ kl ' ks ~ + + W2 77 Kl ' n K2 77,43 ' 7r * 4 

h f 

+ - J d (4) K 7T Kl • 7T K2 7T K3 ■ 7T K4 (3) 

V r 

-— J dK7T K - tt_„ + 0(g 6 ) 
where the inverse propagator, given by 



to 2 

A (/) (k, u) = k 2 + — + f(u) + gh 



depends on the damping term f(co), and the measure term d^n is 

d (2n) /t = {2Tl) d+1 5 {d+1) {K l + ■■■ + K 2n ) d«i • • • dK 2 



^2n 
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In (3) the term proportional to the extended phase-space volume 

S d 2VL r n dco r d d k 



V = 



J-n 2vr J k <i (2vr) d J 



(2n) d + 1 d J-n 2tt J k <i (2tt) £ 
comes from the perturbative evaluation of the Jacobian 

n 1 



2 

X.T 



(4) 



(5) 



In (4) Q is an arbitrary frequency cut-off which does not enter the final result, 
and S d = 2it d / 2 /T(d/2). 

Let us denote the vertices appearing (in the same order) in (3) as 




1 2 

+ 

4 3 




1 2 

+ 

4 3 




1 2 

+ 

4 3 




where the cut, double cut and triangle indicate respectively a k, uj and f(cu) 
insertion, and the dotted line stays for the first perturbative term of the Ja- 
cobian. Let then S — S + S 1: with 



So = Y J dK A (/) (k, u) 7T K ■ 7T_ 



2g 

We shall use the momentum-shell method to derive recursion relations a la 
Wilson [13,40], integrating out the 7r K with e~ l < k < 1, I ~ and all uj, and 
rescaling 

k — > e~ l k, uj — > e~ zl uj, it — > (it 

hereby introducing a dynamic scaling exponent z and a wave function renor- 
malization (. We average with respect to the following propagator: 



(<<>) 



A U) (k,uj) 



, for e 1 < k < 1, e ' < k' < 1 

, for k < e~ l , k' < e~ l 
, otherwise 



The first perturbative term of the Jacobian can be written as 



~ / dKlT K -lT^ K = 

2g Jk<i 2g 



/ dn + I dn dn ir K ■ tt 

Jk<e~ l Je- l <k<l J Jk<l 



(6) 
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and we will indicate it pictorially as 

< > 
• = • + • 



After integration on the momentum shell the first term represents the modified 
first-loop Jacobian contribution, while the second has to be absorbed in the 
mass renormalization and cancels the contribution from the second vertex in 
Eq. (3) (see below), thus maintaining the theory massless at h — 0. 

The usual diagrams renormalizing the k 2 and u 2 terms both give a ^gll(f) 
contribution, with 



S d 



f. 



+°° du 



(7) 



(/) (27^7-00 A (/) (l,^) 
In order to renormalize / and h we have to consider the following contributions: 



- [ d/ci / dn 2 

l Jkx<e- 1 Je- l <k 2 <l 



+ 



(uj 2 /c) 2 | f(ui + u 2 ) ] 1 (N + l)gh 1 



A( f )(k 2 ,u 2 ) A {f) (k 2 ,u 2 ) A {f) (k 2 ,u 2 ) 2A {f) (k 2 ,u 2 



(8) 



2 Jk 1 <e- 1 



- dK 2 + — — ghll {f) 



i<k 2 <\ A(j)(k 2 ,uj 2 



The symmetry of (2) requires that h scales as h — > (h. Since h plays in (3) 
the role of a mass, from one-loop mass renormalization we get that the field 
rescaling factor is 



C c± exp 



N-l T 

l[d + z z—gi(f) 



Defining f(uu) = uj" /c 2 + f(u) we finally get the RG equations: 



dg 
dl 



with 



2 -d-z + (N-2)gI {f) g 



2 - gl {f) - zuJ—j f(u) + gC (f) (uj) 



(9) 
(10) 



'11) 



10 



Equation (9) results as usual from field rescaling and renormalization of the 
k 2 term. The term zu-4- in Eq. (10) describes the frequency rescaling and 
—gl(f) + gC(f) is related to the contributions in Eq. (8). 

Let us now consider a damping term of the form f{uj) = fi(\u\), with f\ a 
generic different iable function behaving as fi(iv) ~ u a for uj — > +oo. One 
observes that for a < 3 and f" integrable in a neighborhood of u = the term 
C(/)(cj) is twice continuously differentiate, so that for uj ~ it is possible to 
expand 

C(/)(cj) = (y\l(f) + 2i^(/)^ ^ 2 + higher order terms (12) 



with 

K - 1 ^ 

{/) 2 (270^ 



/;(») + r i + #^ ? ? \ 

Jo I + UJ Z / c z + fi(u)) 



(13) 



In order to prove (12,13) it is enough to rewrite (11) as 



C(f)H -j2^Jo l + c' 2 /c 2 + A(^) ^ 

Sd ^ h{uj' + uj) + h{uj'-uj)-2f l {uj') , 
+ (2n) d + 1 L 1 + uj' 2 /c 2 + h{u>) U 

and to differentiate the expression twice with respect to uj. C"^(u) is contin- 
uous but not differentiate in uj = 0. 

The content of Eq. (12) is that the non analyticity of the damping term in 
uj = is not reproduced under renormalization. This implies that the f(u>) 
term renormalizes only according to its bare dimensions and to the contribu- 
tion of the field rescaling factor. On the other hand, f(uj) contributes to the 
renormalization of the uj 2 term. 

According to Eq. (13), K(j) is essentially a measure of the non linearity of / 
and is not strongly dependent on its explicit form. 

For f{uj) = r y\uj\ a , with 1 < a < 2 and 7 a damping coefficient, equation (10), 
together with (12,13), gives 



^ = [2-az-0j a (c7)]7 ( 14 ) 



, 2(l-z)- + 2gK a (c n ) (15) 



with 
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-w " (2n) d+1 J-oc l + u 2 /c 2 + -f\cu\ a v 7 

S d 7 /•+°° w a ~ 2 dcu 

^^^(^r^-^io i + ^/c 2 +7^ (17) 

In Appendix A the limit for a — > 1 of these expression is discussed (see 
(A.3.A.4)). 



3 Fixed point analysis of the RG equations 



3.1 One-loop analysis 



We shall now study the RG equations obtained in the previous section for the 
zero-temperature critical point of the quantum non-linear sigma model with 
a damping term f(ou) = 7|w| a , described by the action (2). For the sake of 
clarity we rewrite here the RG equations (9,14,15): 



= [ 2 _ d - z + g(N - 2)I a (c, 7 )]<7 (18) 
^ = [2-OLZ-gI a (c,i)Yt (19) 
~ = 2{l-z)^ + 2gK a {c, 1 ) (20) 

where the integrals I a and K a were defined in (16,17). The dynamic scaling 
exponent z was introduced when choosing to rescale frequencies as uj — > e~ zl cu, 
I being the parameter of the renormalization group. 

We look for fixed points of the RG transformation (18-20). One immediately 
finds two fixed points, characterized by 

FP1: 0*1 = — - 1 , 7* = 0, c* = const, z x = 1, for d > 1; 

(N - 2J7TC* 

ad -2a + 2 2N - 2 - d 

FP2: (7* 2 — 77TT 7T{ , n ,- / 7; c *7* — M *> ^2 — ttt ^ — , 

l(N — 2)a + lj/ a (c*, 7*) (iV-2)a+l 

for d > 2 - -; 

a 

where u* is the solution of 

'•<«> - (21) 
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with I a (u) = c 1 I a (c,uc a ) = I a (l,u), K a (u) = cK a (c,uc a ) = K a (l,u). 

Strictly speaking, FP1 and FP2 are two lines of fixed points, since in both 
cases c* is arbitrary. In the case of FP2 either c* or 7* can be arbitrarily chosen, 
and the remaining one is then determined through the constraint c"7* = u*. 

Concerning Eq. (21), note that, letting fi = 2N ^ d Sf® { ~^__ 2 - )a ■ 

(1) a unique solution u* exists for \i > 0; 

(2) u* — > +00 for /i — > 0; 

(3) no solution exists for \x < 0. 

Together with d > 2 — 2/a, this implies d < 2N — 3 — (N — 2)a. These are all 
consequences of the following properties: I a {u) is a monotonically decreasing 
function of u, tending to a finite value for u — > and vanishing for u — > +00; 
K a (u) is an increasing function, vanishing for u — > and unbounded for 
■u — > +00. 

The analogy with the statistical mechanical model shows that the perturbative 
expansion is in the parameter e = d + — 2, where is the value of the 
dynamic scaling exponent obtained from purely dimensional considerations. 
For the 7 = case Lorentz invariance imposes z^ = 1, while in the presence 
of a non-zero damping term ~ \u\ a with a not too close to 2 we expect the 
latter to be more relevant than uj 2 , and z^ — 2/a. Letting then 5z = z — z^ 
be the one-loop perturbative correction to the bare dynamic exponent z^ and 
rj the field anomalous dimension, one finds in the two cases : 

FP1: *<°> = 1, ^ = 0, Cl = d-1, 9*,i= (N \ c: Vi = 

FP2: 4 0) = -, Sz2 = - % e 2 = d + --2, 

a (N — 2)a + 1 a 

a €2 (X62 

9 *' 2 = (iV-2)a + l/ Q (c*,7»)' V2 = (N-2)a + V 

with c°7* = the solution of I a (u) = -— — — — -K a {u). 

[1\ — Z)a + 1 Z2 — 1 

The point FP2 exists for 62 > 0. For €2 < one gets the unphysical situation 
< 0. The additional condition d < 2N — 3 — (N — 2)a amounts to Z2 = 
a ~ (Af_2 2 )a+i > ^' wn i cn is the condition of stability to order €2 of FP2 with 
respect to FP1. More precisely, at Z2 — 1 FP2 merges into FP1 (see below). 
For £2-^0 one has c"7* — > +00; keeping 7* fixed this gives c* — > +00, as 
expected, since for FP2 one has \jc\ = at 0-th order in €2- For a > 1, 
(16) shows that ^(+00,7) is finite, so that g*^ — 0(e 2 )- The same is true in 
the a = 1 case, but some care must be exerted. A direct computation using 
(A. 3, A. 4) leads to ambiguous results, since in this case the frequency-wave- 
vector integration shell becomes singular: this is discussed more at length in 
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Appendix A. We just observe here that for a — > 1, e 2 — > implies d — > 
and one has J Q ,(+oo,7) ~ One can send £2 — * either with fixed a ^ 1, 
or with fixed d ^ 0, and in both cases the loop integral / a ( c *>7*) is finite at 
0(e 2 ). 

In order to compute the critical exponent v we linearize the RG equations (18- 
20) and find the eigenvalues of the linearized transformation (see Appendix A 
for more detail): 

FP1: ^-= ( jf = e 1 , 0, u^ = 2-a-jp-^ for d > 1; 

FP2; 1=^ = 62, 0, ^ 2 A = 2--+ ^ 6 2 

z/ 2 a (A — 2)a; + 1 

for d > 2 - 2/a, 2A - 3 - d - (A - 2) a > 0. 

This shows that the point FP2 always has one attractive direction in its region 
of existence, while > (< 0) depends on 2A - 3 - d - (A - 2) a < (> 0). 



We can now assess the stability of the two fixed points. In the (d, a) plane we 
distinguish the regions (see Fig. 1 for the A = 3 Heisenberg case): 

FP1: ei > 0, 2N-3-d- (N - 2)a < 0; 

FP2: e 2 > 0, 2A - 3 - d - (A - 2) a > 0, i.e. z 2 > I. 

In region FP1 there exists only the corresponding fixed point, which is stable. 
In region FP2 there appears a second, stable fixed point, while the point FP1 
becomes unstable. The crossing between the regions FP1 and FP2 takes place 
for 2A — 3 — d — (A — 2)a = 0. Notice that on the border line one gets 
Vi = V2 = 2 — a and z 2 — z± — 1, so that when moving from the region FP2 
to FP1 the critical exponents i] and z vary continuously. 

Finally we report in Fig. 2 a numerical computation of the structure of the 
RG flow for typical values of the parameters. To this purpose it is convenient 
to introduce the rescaled variables 7 = c a 7, g = eg, so that (18-20) take the 
form 

d+~g[(N-2)I a (V-K a m}~9 (22) 
a - g(I a (>y) + aK a tf))]7 (23) 

where (22) and (23) are decoupled from (24) and do not contain the dynamic 
exponent z. This means that in the new variables the RG trajectories project 



&g_ 
dl 



{1- 



^ = [2 
dl [ 



_d 1_ 

dZc* 



2[1- 
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Fig. 1. The shaded areas FP1 and FP2 represent the regions of stability of the two 
fixed points in the (d, a) plane in the N = 3 case. 



on the (g, 7) plane and the two lines of fixed points project on two correspond- 
ing fixed points. This fact was not a priori obvious and greatly facilitates the 
analytical and qualitative study of the RG trajectories. 

We notice the interesting fact that there exists a range of values a* < a < 2 
for which one obtains a z = 1 dynamic scaling exponent in the proximity of 
the quantum critical point, even in the presence of a damping term. In other 
words, there can be "soft" damping terms which do not spoil the spin-wave 
behavior of the fixed point. Within the e-expansion for N = 3 a* is given by 
a* = 3 — d (cf. Fig. 1). We shall argue about the extension of these results to 
the d = 2 case in the next subsection. 
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0.2 0.4 0.6 0.8 1 

S d /(2K)^ + g 

Fig. 2. The RG flow projected on the compactified (g, 7) plane (d = 1.2, a = 1, 
TV = 3). 

Beyond the one-loop analysis 



The above analysis of the stability of FP1 vs. FP2 bears a strong similarity 
with the study of the 4 theory with long range interactions [36,37]. For small 
e 1 and e 2 one has a ~ 2 and d ~ 1 and therefore 77 ~ 2 — a, coherently with 
[36]. The two fixed points studied in [37] correspond to our FP1 and FP2; 
in our approach they are selected by the value of z^°\ The continuity of r\ 
for a — * 2 appears as a consequence of their relative stability, as in [37]. The 
computations here are considerably simpler, since first-loop calculations are 
sufficient to establish the general picture. Another difference is in the presence 
of the dynamic exponent z, which was absent in the classical works [36,37]. 

It is worth stressing the fact, noticed after Eq. (12), that the non-analyticity 
of the term \u>\ a is not reproduced under renormalization, at least at first 
loop. As a consequence 7 does not get singular corrections by itself, but only 
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scales according to its bare dimension and to the anomalous field dimensionF]. 
The indication coming from [36-38] is that this will be true also at higher 
orders. This would imply that, as an exact result, 7 scales in the critical 
region according to the exponent 

x 7 = 2 — az — rj (25) 

irrespectively from the approximation. The value z 2 can then be obtained 
directly from the condition x 7 = 0, giving 

2-7] 
Z2 = 

a 

which is in agreement with our one-loop analysis and is expected to remain 
valid at higher loop orders as well. 

This exact relation can be used to settle the problem of the crossover from 
FP1 to FP2. This crossover takes place when z 2 = Z\ — 1, that is, using (26), 
when a = a* = 2 — 77(0*), where 77(a) is the anomalous field dimension for a 
given a. Continuity of 77 through the crossover implies ??(«*) = 771 = 772 with 
rji being the critical exponent of the sigma model with zero damping in d+ 1. 
For two dimensional systems, like the cuprates, d+1 = 3. Numerical estimates 
for the critical indices of the 0(3) symmetric 4 field theory obtained from 
summed perturbation series at fixed dimension 3 (see [41] and Table 25.4 of 
[27]) give 77 ~ 0.033 ± 0.004 and a* ~ 1.966 ± 0.004. In this respect the e- 
expansion greatly overestimates 77 and underestimates a* in d = 2 (see Fig. 1). 
This suggests that, as a consequence of this small value of rj, since the relation 
(26) is probably an exact relation, for a near one the physically relevant line 
of fixed point would be FP2. 

It appears therefore that the z = 1 behavior observed in the cuprates can be 
accounted for in this framework only in the hypothesis that some mechanism 
like the reduction of the spectral weight of low-frequency quasi-particle is so 
strong to effectively produce a \u\ a damping term with a very close to 2. 



4 The ordered phase: a scaling approach 

In the ordered phase the effect of holes in the small-doping limit of the cuprates 
acts as a perturbation of the T = antiferromagnetically ordered ground 
state. The description of the interplay of spin and fermion degrees of freedom 

5 More precisely, this discussion and Eq. (25) below refer to the coupling 7/5 of 
the dissipative term in Eq. (2). In the critical region g = g* and the behavior of 7 
coincides with the behavior of 7/5. 
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(26) 



when approaching the QCP from the ordered side presents major difficulties. 
The simple introduction in the effective sigma model of a ^\oj\ a term is not 
satisfactory, since this term would be relevant also in the proximity of the 
ordered critical point g = 0, and would completely destroy the spin-wave 
picture of the ordered ground state. 

In the ordered phase the appearance of a non-zero staggered magnetization 
No is accompanied by the presence of two quasi-particle bands, separated by 
an energy A (for a mean field analysis see e.g. Ref. [28]). As a consequence, 
one would expect the decay of spin-wave excitations in particle-hole pairs to 
be disfavored for |a>| < A, and an effective damping be felt only for |a>| > A. 
Moreover, in an AFI-PM transition A should vanish at the QCP as well as Nq 
does. Were this not the case, since A = in the disordered phase one could 
hardly match the critical behavior on the two sides of the transition. We are not 
in a condition to directly link the scale A to the staggered magnetization N 
since the precise relation, if any, should emerge from the unknown microscopic 
dynamics of the quasiparticlesQ. In this Section we try to get some insight in 
this difficult problem using simple scaling arguments. We will conjecture for 
A a scaling law with a critical exponent y, computable in principle from the 
full theory including also the fermionic degrees of freedom, and we shall try 
to assess the value of y in a self-consistent way, once the value of a for the 
damping in the disordered phase is given. We emphasize that this procedure 
is not as reliable as the RG computation exposed in the preceding Sections 
and has an heuristic character. 

A model damping term which could take into account the appearance of the 
scale A is 



|6 

(u 2 + A 2 ) : 



2 





\l ,\ b 




r 
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u 


U}\ , 


M 



/(") = 7 , „ ^{ ' ' 7 ' (27) 



with l<a<2<6<3, and r = 7A a_b . We shall later assume that 
A vanishes at the QCP, so that the crossover between the two regimes is 
continuously shifted while approaching the transition. Here we consider A as 
a free parameter at the same level as r and 7. Assuming that changing the 
scale, f(ui) approximately maintains the form (27) with varying values of r 
and 7 one finds from (10) the characteristic exponents 

x 7 = 2 — 77 — az, x r = 2 — 77 — bz, x^ = z. (28) 



Our phenomenological assumptions suggest that damping be relevant only 

6 Notice for instance that in a large-U Hubbard model, the size of the charge gap in 
the insulating phase is not only related to magnetism, but is instead mainly induced 
by the strong correlation 
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at high frequencies. In this regime the spin-wave velocity c sw goes as c sw ~ 
s i-22 Q Q n o^gj- hand, at any given scale the relevant cj-modes are 
characterized by the condition uo 2 = 0(1). So, damping is relevant up to the 
scale s where A becomes of order 1. Beyond this scale the gap for spin- wave 
decay is visible and f(oo) is irrelevant. If A at the "microscopic" scale s = 1 had 
a finite value, eventually f{u>) would be always irrelevant. The simplest way 
to take into account the physical condition that A should vanish approaching 
the transition is to assume that 

&\ s =i~(9*-9o) v (29) 

where go = g\ s =i an d 9o — ► 9*- This condition encodes in a simplified, phe- 
nomenological way the complicated interaction between damping and anti- 
ferromagnetic order, mediated in the full theory by the fermionic degrees of 
freedom. 

Taking into account (29) and the scaling exponents (28) we get 

s~(g*-g r y/z \ go^g* (30) 

Approximately at this scale damping becomes irrelevant and the critical be- 
havior starts being characterized by z — 1. This gives c sw vanishing as 

Csw ~ s l ~ Z2 ~ (<7* - g ) y{Z2 - 1)/z \ g - g* (31) 

For a not too close to two (specifically for a < a*) by approaching the quan- 
tum transition from the disordered phase one observes, as already discussed, 
a dissipative behavior dominated by the fixed point FP2, with the scaling law 
uj ~ k Z2 . Deep into the broken symmetry phase one instead has a propagat- 
ing behavior, uo ~ c sw k. This should match the critical behavior at the scale 
£ ~ (go — g*) u , with v the corresponding critical index. Continuity of the crit- 
ical behavior in a neighborhood of the critical point imposes then c sw ~ £ 1-22 , 



Let us briefly comment about the relation of the physical spin-wave velocity c sw 
with the running coupling constant c and the dynamic exponent z. In the neighbor- 
hood of a critical point, z is determined by the condition that c does not renormal- 
ize. This corresponds to attributing all the renormalization of the uj 2 /c 2 term to the 
anomalous frequency scaling, lo ~ k z . The running coupling constant c differs from 
the physical spin- wave velocity c sw , which is measured at any scale s using a fixed 
system of physical units: one has c sw ~ s 1_z c. So, although c does not renormalize 
by construction, when z > 1 the physical velocity c sw vanishes at the transition. In 
the case of the undamped sigma model the value z = 1 is obtained as a consequence 
of Lorentz invariance, so that there is no difference between c and c sw and the latter 
remains finite at the QCP. However, the damping term breaks Lorentz invariance 
and can drive c sw to zero at the QCP. 



19 



i.e. 

c sw ~ \g* - 0or (z2-1 \ 9o -> 
This is compatible with (31) only for 

y = ^2 (32) 

In terms of £ this condition on the critical behavior of the energy separation 
A gives 5 ~ i~ Z2 - 



5 Conclusions 

It is not unreasonable to expect that the damped sigma model we considered 
contains some of the physics of the quantum critical point of the underdoped 
layered cuprates. For this reason we performed a RG analysis directly on 
this model, trying to maintain sufficient generality in the form of the effective 
damping term, in particular introducing an exponent a > 1 that could account 
for the observed pseudogap behavior. In the broken symmetry phase we tried 
to get some insight on the interplay between damping and AF ordering using 
simple scaling arguments. 

In both the ordered and the disordered phases, assuming continuity of the hy- 
drodynamic behaviors at the quantum transition, we obtained that for a < a* 
(where a*(d = 2) is expected to be very close to 2) it is not possible to account 
for the observed z — 1 scaling law. It was observed in Ref. [28] that the z — 1 
scaling can be recovered at higher temperatures, when typical frequencies 
would become larger than the T = value of the damping term. However, ex- 
perimentally it is the z — 1 behavior which is observed at lower temperatures, 
while a relaxational z = 2 behavior is observed at higher T [23]. It was argued 
[42,23] that this could be due to the neglected T-dependence of the damping 
term 7. We propose in this paper that the same effect can be accounted for by 
assuming that the decay of the quasi-particle spectral weight for low frequen- 
cies produces a "soft" damping term ~ \uo\ a with a > a*. The z — 1 scaling 
would then be a consequence of an almost linear decay of the quasi-particle 
spectral weight. We notice that this reduced weight is fully consistent with 
the observed occurrence of a pseudogap below a crossover temperature T* 
in the underdoped region of the phase diagram for the cuprates. The origin 
of this pseudogap might arise from extrinsic mechanisms like local Cooper- 
pair formation and/or stripe fluctuations. Alternatively a pseudogap could be 
the very consequence of the intrinsecally non Fermi liquid nature of the sys- 
tem [34] or of strong critical fluctuations occurring near the QCP. In this last 
regard singular AF fluctuations have been shown within a spin-fermion model 
to strongly modify the quasiparticles properties [43] in the proximity of the 
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"hot" points of the Fermi surface separated by the AF wavevector. Since the 
particle-hole excitations around these points are responsible for the low-energy 
spin-wave decay, this naturally affects the damping of spin excitations, in a 
complicated interplay which deserves further investigation. We notice that the 
fluctuation-induced spectral weight reduction for the fermionic excitations is 
a definite signature of a non-Fermi liquid metallic state. This is in agreement 
with general arguments put forward by Anderson (see e.g., Ref. [44]) on the 
difficulty in continuously connecting a charge-gapped insulating phase with a 
normal Fermi-liquid phase. We remark, however, that Anderson's arguments 
extend the presence of a non-Fermi-liquid phase all over the metallic state as 
a consequence of low dimensionality. On the other hand, within the QCP sce- 
nario, the violation of the Fermi-liquid properties is a consequence of critical 
fluctuations, which are only present around the critical point. As soon as one 
moves away from this point, the fluctuations loose their critical character thus 
providing simple perturbative corrections to the Fermi liquid state. 

It is worth noting that the analysis carried out in this paper is not only re- 
lated to the superconducting cuprates, but it has connections with other phys- 
ically interesting problems, like the spin chains with long range interaction and 
macroscopic quantum tunnelling. By considering the model (2) in d — 0, it 
is seen that the frequency axis is left as the only relevant "direction" (cf. the 
technical considerations presented in Appendix A). As a first consequence, all 
the critical behavior should be reformulated in terms of frequency by suit- 
ably rescaling of the critical indices by \jz. The frequency axis becomes then 
equivalent to a single space direction, via the dynamic critical index z, thus 
establishing a connection with onedimensional classical models. Fourier trans- 
forming back to (imaginary) times the damping term in (2) gives rise to a long 
range interaction of the type \t — t'\~ 1 ~ a . Our results can be directly extended 
to the case < a < 1. This leads to the comparison with the onedimensional 
(classical) iV-component spin model with long-range interaction considered 
in Ref. [45]. One can recognize that, once the spatial term is dropped in (2) 
the connection is established through the coupling combinations 7/5 and gc 2 . 
This allows to rewrite the RG equation for g/7 in the same form as Eq. (7) of 
Ref. [45], with g/j (to be identified with T) having a fixed point (FP2 in our 
model) and (l/gc 2 )* = 0. 

Whitin the d = case, N = 2 deserves special attention because of the connec- 
tion with macroscopic quantum tunnelling [46,47]. This connection is realized 
by observing that the 0(2) non-linear a model can be rewritten in terms of an 
angular degree of freedom and that the addition of the magnetic field produces 
a cos-like interaction (sine-Gordon model). In d — the model then describes 
a single degree of freedom with a kinetic (u 2 ) term in a multivalley potential. 
In the presence of dissipation, a linear-in-frequency term can appear, leading 
again to long-range interactions on the (imaginary) time axis. Models of this 
kind were considered by Chakravarty [48] in the case of a double well poten- 
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tial and successively by Schmid [49] to describe a dissipative quantum particle 
in a onedimensional periodic potential. In particular, the only difference be- 
tween our Eq. (2) at N = 2 and this latter model is the non-periodicity of the 
dissipative term of Ref. [49]. 

In the model of Schmid, for vanishing h (g = in his notation), a line of fixed 
points for 77 (corresponding to our j/g) was found, with a critical point r\ c 
separating a regime with relevant h from one with irrelevant h. The long-range 
term introduced in [49] is quadratic in the fields, so that the model for h — is 
free and r\ (i.e. j/g) does not renormalizeQ. A (dual) fixed point was also found 
for very large h, where the model of Ref. [49] was stated to be equivalent to the 
one of Ref. [48]. We observe that, differently from the model of Ref. [49] our 
long-range term mantains the periodicity in the angular variable and provides 
an interaction. This renormalizes the coupling, thus allowing j/g to flow under 
renormalization even at h = 0. One could conjecture that a new fixed point 
arises at finite h. An adequate analysis would require a RG study at finite h 
going beyond the e-expansion and taking into account instantonic corrections. 
This is an interesting open problem which deserves further investigation. 
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A Computation of the critical indices 



For the computation of the eigenvalues u R , u A of the linearization of the RG 
transformation (18-20) in the point FP2 it is convenient to rewrite it in the 
form (22-24). The secular equation is found to depend on 



?a = u-^-\ogI a (u) 



(A.l) 



with u — 7 = c°7 and I a and K a given after Eq. (21). For e = d — 2 + 2/a — > 
one has u* — > +oo, and in this limit it is easy to derive from (16,17), for a > 1: 

Ia[U) ~ (2n) d +i asm ^ ' KaW ~ (2tt)^ 2sin^ U 
which substituted in (A.l) gives 

T a ~ --, K a ~ - (A.2) 
Solving (21) asymptotically for -u* — > +oo (e — > 0) and fixed 7* gives 



( 4[(A^-2)q + l](2-ct) | Q 1 J_ 
\ aP(a - 1) J C *^7^ 



Using (A.2) one easily finds o;^, cj^ at first order in e. 
For a — > 1 the expressions (16,17) reduce to 




- .... ...... V / ''-( f ") 2 



(27T) d+1 2 



v/I -^arctan^ T -,0<c 7 <2 

1 og ^^ , C7>2 
V( c 7) 2 -4 6 c 7 - v /( c 7) 2 "4 ' ' 

(A.4) 



The same expressions can be obtained by a direct evaluation of (7, 13) for a = 
1. In particular the original expression (7) shows that I\ is in fact an analytic 
function of c and 7. However, some problems arise when using expressions 
(A.3) and (A.4) to compute critical indices directly for a = 1 instead of 
considering the a — > 1 limit of the general expressions given in Section 3.1. 
For a — 1 in fact one expands around the critical dimension d c = and the 
factor Sd coming from angular integration vanishes. If the same computation 
is performed using field-theoretical methods, instead of Wilson's integration 
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on the momentum shell, no such problems arise: e.g. one has 



/ 



k 2 + u 2 /c 2 + ^uj- 



dcu d d k 



2 7T/2 2^ 2 r (|5f) r (^y) 

2-asinf r(f) r(l-f) 



ae d—1 



(A.5) 



which correctly behaves as e 1 for e — > 0, independently of the value of a. 
Notice that the r _1 (d/2) factor coming from Sd is cancelled by the factor 
(sin(d7r/2)) _1 and the e -1 divergence arises from r(ae/2(2 — a)). The problems 
with the Wilson method are resolved if one makes a more symmetric choice of 
the integration shell, e.g. integrating over all the 7r fc w with e~ l < \Jk 2 + oo 2 < 1 
(which implies introducing a frequency cutoff). With this choice I a is finite for 
d — > 0. We observed that I a , K a computed with either choice of cutoffs have 
the same c, 7 dependence, the only effect of the choice of the symmetric shell 
being that the integrals do not vanish for d — > 0. In particular, the expressions 
for uf, 0J2 given in the text are valid for 1 < a < 2: as a matter of fact, 
they are computed from the values (A. 2), which do not contain the angular 
integration factor Sd. 
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